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Abstract 

An acoustic cloak envelopes an object so that sound incident from all directions passes 
through and around the cloak as though the object were not present. A theory of acoustic 
cloaking is developed using the transformation or change-of- variables method for mapping 
the cloaked region to a point with vanishing scattering strength. We show that the 
acoustical parameters in the cloak must be anisotropic: either the mass density or the 
mechanical stiffness or both. If the stiffness is isotropic, corresponding to a fluid with a 
single bulk modulus, then the inertial density must be infinite at the inner surface of the 
cloak. This requires an infinitely massive cloak. We show that perfect cloaking can be 
achieved with finite mass through the use of anisotropic stiffness. The generic class of 
anisotropic material required is known as a pentamode material. If the transformation 
deformation gradient is symmetric then the pentamode material parameters are explicit, 
otherwise its properties depend on a stress like tensor which satisfies a static equilibrium 
equation. For a given transformation mapping the material composition of the cloak is 
not uniquely defined, but the phase and wave speeds of the pseudo-acoustic waves in the 
cloak are unique. Examples are given from 2D and 3D. 

1 Introduction 

The observation that the electromagnetic equations remain invariant under spatial transfor- 
mations is not new. Ward and Pendry [1] used it for numerical purposes, but the result 
was known to Post [2] who discusses it in his book, and it was probably known far earlier. 
The recent interest in passive cloaking and invisibility is due to the fundamental result of 
Greenleaf, Lassas an Uhlmann [3l 1^ that singular transformations could lead to cloaking for 
conductivity. Not long after this important discovery, Pendry et al. [5] and Leonhardt [6] 
made the key observation that singular transformations could be used to achieve cloaking of 
electromagnetic waves. These results and others have generated significant interest in the 
possibility of passive acoustic cloaking. 

Acoustic cloaking is considered here in the context of the so-called transformation or 
change-of-variables method. The transformation deforms a region in such a way that the 
mapping is one-to-one everywhere except at a single point, which is mapped into the cloak 
inner boundary, see Figure [H The acoustic problem is for the infinitesimal pressure p(x, t) 
which satisfies the scalar wave equation in the surrounding fiuid, 

V^p -p = 0. (1) 

The basic idea is to alter the cloak's acoustic properties (density, modulus) so that the 
modified wave equation in uj mimics the exterior equation ([T]) in the entire region Q. This is 
achieved if the spatial mapping of the simply connected region Q to the multiply connected 
cloak UJ has the property that the modified equation in lo when expressed in 0, coordinates 
has exactly the form of ([T]) at every point in $7. 
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Figure 1: The undeformed simply connected region is transformed by the mapping x iiito 
the multiply connected cloak uj. Essentially, a single point O is transformed into a hole (the 
invisible region) surrounded by the cloak to. The outer boundary is coincident with 
dQ+{= dQ) and the inner boundary dco- is the image of the point O. Apart from O and 
duj- the mapping is everywhere one-to-one and differentiable. 

The objective here is to answer the question: What type of material is required to realize 
these unusual properties that make an acoustic cloak? While cloaking cannot occur if the bulk 
modulus and density are simultaneously scalar quantities (see below), it is possible to obtain 
acoustical cloaks by assuming the mass density is anisotropic [71 El [9] . A tensorial density 
is not ruled out on fundamental grounds [TO] and in fact there is a strong physical basis for 
anisotropic inertia. For instance, Schoenberg and Sen [11] showed that the inertia tensor in a 
medium comprising alternating fluid constituents is transversely isotropic with elements {p) 
in the direction normal to the layering, and in the transverse direction, where (•) 

is the spatial average. Anisotropic effective density can arise from other microstructures, as 
discussed by Mei et al. [12] and by Torrent and Sanchez-Dehesa [13]. The general context 
for anisotropic inertia is the Willis equations of elastodynamics jl4j which Milton et al. [lOj 
showed are the natural counterparts of the EM equations that remain invariant under spatial 
transformation. Acoustic cloaking has been demonstrated, theoretically at least, in both 2D 
and 3D: a spherically symmetric cloak was discussed by Chen and Chan [8] and by Cummer 
et al. [9], while Cummer and Schurig [7] described a 2D cylindrically symmetric acoustic 
cloak. These papers use a linear transformation based on prior EM results in 2D [15| . 

Cloaks based on anisotropic density in combination with the inviscid acoustic pressure 
constitutive relation (bulk modulus) will be called Inertial Cloaks (IC). The fundamental 
mathematical identity behind the IC is the observation of Greenleaf et al. |16| that the scalar 
wave equation is mapped into the following form in the deformed cloak region 



Here g = (gij) is the Riemannian metric with \g\ = det{gij), (5*-') = {gij)~^. The reader 
familiar with differential geometry will recognize the first term in eq. (|2]) as the Laplacian in 
curvilinear coordinates. Comparison of the transformed wave equation ([2]) with the IC wave 
equation provides explicit expressions for the IC density tensor and the bulk modulus [T7] . 

We will derive an identity equivalent to ([2]) in Section [2] using an alternative formulation 
adapted from the theory of finite elasticity. A close examination of the anisotropic density of 
the IC shows that its volumetric integral, the total mass, must be infinite for perfect cloaking. 
This raises grave questions about the utility of the IC. The remainder of the paper provides 
a solution to this quandary. The main result is that the IC is a special case of a more general 
class of acoustic cloaks, defined by anisotropic inertia combined with anisotropic stiffness. 
The latter is obtained through the use of pentamode materials (PM) [18]. In the same way 
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that an ideal acoustic fluid can be defined as the limit of an isotropic elastic solid as the shear 
modulus tends to zero, there is a class of limiting anisotropic solids with five (hence penta) 
easy modes of deformation analogous to shear, and one non-trivial mode of stress and strain. 
The general cloak comprising PM and IC is called the PM-IC model. The additional degrees 
of freedom provided by the PM-IC allow us to avoid the infinite mass dilemma of the IC. 

We begin in Section [2] with a new derivation of the inertial cloak (IC) model, and a 
discussion of the infinite mass dilemma. Pentamode materials are introduced in Section [3] 
where it is shown that they display simple wave properties, such as an ellipsoidal slowness 
surface. The intimate connection between PM and acoustic cloaking follows from Theorem [1] 
in Section m Properties of the generalized PM-IC model for cloaking are developed in Section 
m through the use of an example cloak that can be pure IC or pure PM as a parameter is 
varied. Further examples are given in SectionEJ with a concluding summary of the generalized 
acoustic cloaking theory in Section [6l 

2 The inertial cloak 

The transformation from to u; is described by the point-wise deformation from X G O 
to X = x(X) G uj. In the language of finite elasticity, X describes a particle position in 
the Lagrangian or undeformed configuration, and x is particle location in the Eulerian or 
deformed physical state. The transformation or mapping defined by x is one-to-one and 
invertible except at the single point X = O, see Figure [TJ We use V, Vx and div, Div to 
indicate the gradient and divergence operators in x and X, respectively. The component form 
of div A is dAi / dxi or dAij / dxi when A is a vector and a second order tensor-like quantity, 
respectively. The deformation gradient is defined F = Vxx with inverse F~^ = VX, or 
in component form Fn = dxi/dXj, Ff-^ = dXj/dxi. The Jacobian of the deformation is 
J = detF = |F|, or in terms of volume elements in the two configurations, J = dv/ dV. 
The polar decomposition implies F = VR, where R is proper orthogonal (RR* = R*R = I, 
detR = 1) and the left stretch tensor V GSym"*" is the positive definite solution of = FF*. 
The analysis is as far as possible independent of the spatial dimension d, although applications 
are restricted to d = 2 or d = 3. 

The principal result for the inertial cloak is 

Lemma 1 

V^p = J div ( V^Vp) . (3) 

Proof: The right hand side can be expressed 

J div ( J-^FF*Vp) = J J"^ (F*V) • (F*Vp) + J(F*Vp) • div ( J~^F) . (4) 

Using the chain rule in the form F*V = Vx or V = F~*Vx implies that F* div (F*Vp) = 
Div VxP which is V^p. The proof follows from the identity (see Problems 2.2.1 and 2.2.3 in 

m) 

div (J"^F) = 0. (5) 
2.1 Cloak acoustic parameters 

The connection with acoustics is made by identifying the field variable p in Lemma [T] as the 
acoustic pressure. The cloak comprises an inviscid fluid with bulk modulus K{x.) such that 
the pressure satisfies the standard relation 

p = —Kdivv, (6) 

where v(x, t) is particle velocity. The inertial cloak is defined by the assumption that the 
momentum balance involves a symmetric second order inertia tensor p according to 

pv = -Vp, (7) 
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Although this is a significant departure from classical acoustical theory in assuming an 
anisotropic mass density, it is by no means unprecedented. Based on the analysis of Schoen- 
berg and Sen [11], a spatially varying tensor p could possibly be achieved by small pockets 
of layered fluid separated by massless impermeable membranes. 

Eliminating the velocity between eqs. ^ and ([7]) gives a single equation for the pressure 

K div [p'^Vp] - p = 0, X e w. (8) 
Consider the uniform wave equation in Vt: 

Sj\p-p = {)^ XgQ. (9) 

Using Lemma [J we can express this in the deformed physical description as eq. ([8|) where 
the bulk modulus and inertia tensor are 

K = J, p = JV~2 (;L0) 

For a given deformation F, the identities (jlOp define the unique cloak with spatially varying 
material parameters K and p each defined by the deformation gradient. We note the following 
identity which is independent of F: 

det p = K'^^\ (11) 

Could the cloak possibly have isotropic density? That is, could the cloak be described by 
a standard acoustic fluid with two scalar parameters, density and bulk modulus? The identity 
p = JV~^ means that p = pi can occur only if V is a multiple of the identity, V = for 
some scalar w = w{x). The deformation of 0, into the smaller region u! could certainly be 
accomplished at some but not all points by this deformation, which corresponds to a uniform 
contraction or expansion, with rotation. However, the deformation near the inner surface of 
the cloak cannot be of this form. In fact, the deformation in the neighbourhood of X = O 
must be extremely nonuniform and anisotropic. We will discuss this below when we examine 
a fundamental and severe deficiency of the inertial cloak model. 



2.2 Continuity between the cloak and the acoustic fluid 

Let ds, n and dS", N denote the area element and unit normal to the outer boundary duj^, 
9f2_|_(= duj^), respectively. These are related by the deformation through Nanson's formula 
[l9]NdS'=J"^F*nds. The nature of the cloak requires that the outer surface is identical 
in either description, since both must match with the exterior fluid. We therefore require 
that d s = d S" at every point on the outer surface, or 

N = J-^F*n on dLo+ = dn+, (12) 

and eq. (flO]) then implies 

p-^n = FN on duj+ = dn+. (13) 

Equation (jl3p is a purely kinematic condition. 

The interior of the cloak mimics the wave equation in the exterior fluid. The flnal re- 
quirement that the cloak will be acoustically "invisible" is that the pressure and normal 
velocity match across the outer surface separating the fluid and cloak. These two continuity 
conditions arise from the balance of force (normal traction) per unit area and the constraint 
of particle continuity. The condition for pressure is simply that p is continuous across the 
outer surface, whether one uses the wave equation in physical space, ([8]), or its counterpart 
in the undeformed simply connect region, ([9]). As for the kinematic condition, consider its 
equivalent, the continuity of normal acceleration. This is t'n = n • v in physical space, and 
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using eq. ([7]) it becomes Vn = — n • p ^Vp, which must match with — n • V^* in the fluid. 
Alternatively, eq. ([T3|) and the relation F*V = Vx imply, as expected, 

Vn = -n- p'^Vp = -NF* • Vp = -N • VxP- (14) 

The final term is simply the normal acceleration in the undeformed description. 

In summary, the continuity conditions at the outer surface in the physical description are 

[p]=0, [n-p-^Vp\={) onduj+. (15) 

2.3 Example: a rotationally symmetric inertial cloak 

Consider the inverse deformation 

X = /(r)x, (16) 
where x = x/r and r = |x|. Using F^^ = VX implies 

F = (l//Ol. + (r//)I^, (17) 

where /' = d//dr and the second order tensors are I,. = x x, = I — x x. The bulk 
modulus and mass density in the cloak follow from eq. (jlOp as 



The anisotropic inertia has the form 

p = Prlr + PAj±, (19) 



where the radial and azimuthal principal values and p_\_ can be read off from eq. (jlSp as 

functions of /. 

Introducing the radial and azimuthal phase speeds: = ^J'Kfp'r and c_l = Kj p± , the 
mass density tensor can then be expressed p = K{c~'^lr + cJ^I_l). The quantity Kpr is the 
square of the radial acoustic impedance, z^. = y/Kpr- Equation (fTTj) implies that the identity 
required for cloaking. The three equations (fT8]) for K, pr and p± in terms of / 
can be replaced by the universal relation pT]). i.e., 

Prpt' = (20) 

along with simple expressions for the wave speeds in terms of /: 

Cr = J, = J. (21) 

We will see later that the phase and the wave (group velocity) speeds in the principal direc- 
tions are identical. Note that /' is required to be positive. The original quantities can be 
expressed in terms of the phase speeds as 

p^ = c;^c'^-\ p^ = CrC'^-^, K = CrC'^~\ (22) 

One could, for instance, eliminate / as the fundamental variable defining the cloak in favor 
of c±{r), from which all other quantities can be determined from the differential equation 
relating the speeds: {r/c±) = l/cr- 

We assume the cloak occupies w = {x : < a < |x| < 5} with uniform acoustical 
properties K = 1, p = I in the exterior. The areal matching condition (|13p with n = N = x 
is satisfied by F and p of eqs. (I17p and (I18|) if / is continuous across the boundary, which 
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Table 1: Behaviour of quantities near the inner surface r = a for the scahng / oc as 
^ = r — a I 0. The total radial mass rrir is defined in eq. (j24p . 
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is accomplished by requiring f{b) = b. The pressure and velocity continuity conditions (fT5]) 
become 

H=0, [j^]=0 on r = b. (23) 

Note that the cloak density is isotropic if Cr = c±, which requires that /' = f /r. Thus 
f = with 7 constant, but the outer boundary condition f{b) = b implies 7 = 1, which is 
the trivial undeformed configuration. 

Perfect cloaking requires that / vanish at r = a. It is clear that c± blows up as r | a, as 
does the product Kpr- In order to examine the individual behaviour of K and pr consider 
/ oc (r — a)" near a for a constant and non-negative. No value of a > will keep the radial 
density pr bounded, although the unique choice a = 1/d ensures that the bulk modulus K{a) 
remains finite and non-zero. Note that the azimuthal density p± has a finite limit in 2D for 
power law decay / oc (r — a)", while p± remains finite in 3D iff a < 1, otherwise it blows up. 
Similarly, the radial phase speed scales cis oc (r — a)^ which remains finite for a < 1, 
blowing up otherwise. These results are summarized in Table 1. 

We use a nondimensional measure of the total mass in the cloak: m = (vol(ct;))~^ J dvp. 
The total mass is isotropic for the symmetric deformation and configuration considered here: 
m = ^("T-r + {d — l)m_|_)I where (mr,m_L) = {yoI{lu))~^ f^dv{pr, p±). Assuming for the 
moment that /(a) is nonzero, i.e., a near-cloak [2D], then 



mr={^ t^'i" ^^^^"^ - 2 /: d r- r In fir)] , 2D, ^^^^ 



These forms indicate not only that rtir — > oo as /(a) — > 0, but also the form of the blow-up. 
To leading order, rrir = j^z^ jf^ -|- ... in 2D and = j^z^jf^ -|- ... in 3D. The blow-up 
of nir occurs no matter how / tends to zero. The infinite mass is an unavoidable singularity. 



2.4 A massive problem with inertial cloaking 

Table 1 and the example above illustrate a potentially grievous issue: infinite mass is required 
for perfect cloaking in the IC model. We now show that the problem is not specific to the 
rotationally symmetric cloak but is common to all inertial cloaks. Consider a ball of radius e 
around X = O. Its volume d V =0{e'^) is mapped to a volume with inner surface defined by 
the finite cloak inner boundary duj- and outer surface a distance O(e^) further out, where 
/? > is a local scaling parameter, assumed constant (in terms of the example above and Table 
I, P = 1/a). The mapped current volume is then dv =0(e^) so that J = dv/ dV =0{e^~'^). 
The eigenvalues of V are Ai =0(e^~^), X2,...,d =0(e~^). The bulk modulus and the principal 
values of the density matrix are therefore 

K = 0{e^-% pi = 0(e2-'^-^), P2,...4 = 0{e^-''+'). (25) 

The principal value pi blows up whether d = 2 or d = 3. Furthermore, the total mass 
associated with pi in the mapped volume is mi =0(6^"^^) which blows up in 3D, and a more 
careful analysis for 2D similar to that for the rotationally symmetric case shows mi =0(| In e|). 
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In summary, the inertial cloak theory, while consistent and formally sound, reveals an 
underlying and "massive" problem. We will show how this can be circumvented by using a 
more general cloaking theory which allows for anisotropic stiffness (elasticity) in addition to, 
or instead of, the anisotropic inertia. The anisotropic elastic material required is of a special 
type, called a pentamode material [18], which is introduced next. 



3 Pentamode materials 



We consider Hooke's law in 3D in the form <t 
and the associated 6x6 matrix of moduli are 



Ce, where the 6-vectors of stress and strain. 
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The V2 terms ensure that products and norms are preserved, e.g. CijkiCijki = trC C. 

A pentamode material (PM) is rank one, or in other words, five of the six eigenvalues of 
C vanish [18]. The one remaining positive eigenvalue is therefore 



= tr C = Ci,y = Cu + C22 + C33 + 2(^44 + C55 + Qe). (26) 
Accordingly, the moduli can be defined by the stiffness K and a normalized 6-vector s. 



C = Kss\ 3 

The stress is described by a single scalar, <t = as with a = Ke and e = s*£. Thus, 



(27) 



C = A' S S 



ijkl 



K SijSki, 



V2- 



Si 



(28) 



The pentamode material [T^ is so named because there are five easy ways to deform it, asso- 
ciated with the eigenvectors of the five zero eigenvalues of the elasticity stiffness. Pentamodes 
obviously include isotropic acoustic fluids, for which the only stress-strain eigenmode is a hy- 
drostatic stress, or pure pressure, and the five easy modes are all pure shear. Milton and 
Cherkaev [18] describe how pentamode materials can be realized from specific microstruc- 
tures. 



3.1 Example: An orthotropic PM 

An elastic material with orthotropic symmetry has nine non-zero elements in general: the 
six Cij = Cji, i,j = 1,2,3, plus C44, C55 and Cqq. We set these last three (shear) moduli 
to zero. The stress cr must then be diagonal in the cartesian coordinate system, implying 
S4 = §5 = sq = 0, and therefore 

K = ^{Cu + C22 + C33), S = ^-1/2 (CiY'ei ® ei + Cli^e2 ® 62 + C3'3/'e3 ® e^) , 
with the following relations holding: C12 = c|/^C22^, C23 = ^22^6*33^, C13 = cK'^cl^'^. 
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3.2 Compatibility condition for pentamode materials 

The notation K and S is used to signify the fact that the S tensors are normalized by tr S =3 
and therefore K is given by eq. (j26p . We will not follow this normalization in general, but 
write 

C = irS®S ^ Cijki = KS^,SH. (29) 

In other words, the products in ([29]) are the important physical quantities, not K and S 
individually. The stress in the PM is always proportional to the tensor S and only one strain 
element is significant, S : e. The rank deficiency of the moduli, which is apparent from ()27p 
or ()29p . means that there is no inverse strain-stress relation for the elements of e in terms of 
the elements of a. 

Static equilibrium of a pentamode material under an applied load leads to a constraint on 
the spatial variability of the PM stiffness. Consider an inhomogeneous pentamode material 
with smoothly varying C(x) = i^Tg (x)So(x) So(x). Under an applied static load the strain 
will also be spatially inhomogeneous, but the only part of the strain that is important is the 
component along the PM eigenvector. With no loss in generality we may put e(x) = t(;(x)So 
for some scalar function w. The stress is then cr = gSo where g(x) = wKq tr(So) = "iKw. 
Let S = gSo, then the static equilibrium condition diver = becomes divS = 0. Finally, the 
PM stiffness is C = S where K(x) = KQ/{?,Kwf. 

Lemma 2 The fourth order stiffness of a smoothly varying pentamode material can always 
he expressed C = KS ® S, where K{x) > and S(x) £Sym satisfies the static equilibrium 
condition 

divS = 0. (30) 

This identity also arises in a completely different manner later when we consider transformed 
wave equations. We say that the PM is of canonical form when eq. (I30p applies. The 
decomposition of Lemma [2] is unique up to a multiplicative constant. Thus, if a static load is 
applied to a PM expressed in canonical form, then the stress and strain are <t(x) = cqS and 
e(x) = {co3K)~^S, respectively, for constant cq. 

In summary, stability under static loading places a constraint on the PM moduli, which 
will turn out to be useful when we return to the cloaking problem. The constraint means 
that the moduli can in general be expressed in canonical form. 

3.3 Dynamic equations of motion in a PM 

The equations for small amplitude disturbances in a PM with anisotropic mass density are 

cr = Ktr(Se)S, (31) 
pv = diver. (32) 

These are respectively, the specific form of Hooke's law for a PM and the momentum balance 
incorporating the inertia tensor. In order to make the equations look similar to those for 
an acoustic fluid, we identify the "pseudo-pressure" p with the negative single stress: p = 
—Ktr{Ss). The stress tensor becomes 

cr = —pS, (33) 

and the linear constitutive relation can be written 

p = -KS : Vv. (34) 

Equations (f32]) and (l34|) imply that the pseudo-pressure satisfies the generalized acoustic 
wave equation 

/STS: V(p~Miv(Sp)) -p = 0. (35) 
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This reduces to the acoustic equation ^ with anisotropic inertia and isotropic stiffness when 
S = I. Finally, assuming that the PM is in canonical form, so that S satisfies the equilibrium 
condition ([30]) . we have 

KS : V{p-^SVp) -p = 0. (36) 
3.4 Wave motion in a PM 

The wave properties of pentamode materials are of interest since we will show that they can 
be used to make the acoustic cloak. Consider plane wave solutions for displacement of the 
form u(x, t) = qe*'^('^'-'^~^'*) , for |n| = 1 and constant q, k and v, and uniform PM properties. 
Non trivial solutions of the equations of motion (I3ip - (l32p must satisfy 

{K (Sn) ® (Sn) - pt;2) q = 0. (37) 

The acoustical or Christoffel |21j tensor i^(Sn) (g) (Sn) is rank one and it follows that of the 
three possible solutions for v'^, only one is not zero, the quasi-longitudinal solution 

v'^ = Kn- Sp~^Sn, q = p~^Sn. (38) 

The slowness surface is therefore an ellipsoid. Standard arguments for waves in anisotropic 
solids |21j show that the energy flux velocity (or wave velocity or ray direction) is 

c = v~^KSp~^Sn. (39) 

Note that this is in the direction Sq, and satisfies c-n = v, a well known relation for generally 
anisotropic solids with isotropic density. 

As an example, consider the orthotropic PM with a density tensor of the same symmetry 
and coincident principal axes. Then 

v"^ =c\n\ + clnl + clnl, (40a) 
c =v~^ (cf niei + 0271262 + cfnses) , (40b) 

-1/2 , "1/2 , -1/2 N 

q=Pi ciniei + p2 C2?i2e2 + 0371363, (40c) 
where cf = Cu/pi, c| = C22/P2, c| = C^s/ps, and pi, p2, P3 are the principal inertias. 



4 The general acoustic cloaking theory 

We now show that the inertial cloak (IC) is but a special case of a much more general type 
of acoustic cloak. While the IC depends upon the anisotropic inertia, the general cloaking 
model can have both anisotropic inertia and stiffness. The additional degree of freedom is 
obtained by replacing the pressure field with the scalar stress of a PM. The general cloaking 
model is cahed PM-IC. 

4.1 The fundamental identity 

Lemma 3 Let P G Sym be non-singular and F is the deformation gradient for the mapping 
X ^ x with J = det F, = FF*. Then 

V\p = JP : V(j~^P"^v2Vp), (41) 

iffP satisfies 

divP = 0. (42) 

The proof is given in Appendix A. This clearly generalizes Lemma [U and in the context of 
pentamode materials it implies 
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Theorem 1 The pressure p satisfies a uniform wave equation in Q. Under the transforma- 
tion Q ^ uj with J = detF, = FF*, p satisfies the equation for the pseudo-pressure of a 
pentamode material with stiffness C and anisotropic inertia p: 

V^p-p = inn ^ KS : V(p"^SVp) -p = inuj, (43a) 

where 

K = J, C = KS(g)S, p = JSV-^S, (43b) 

and S satisfies 

div S = 0. (43c) 

Note that the stress tensor S is not uniquely defined, ahhough it must satisfy the equilibrium 
condition (j43cp . The associated density depends only on the left stretch of F, viz., V. The 
inertial cloak corresponds to the special case of S = I, which is a trivial solution of eq. (j43cp . 
The importance of Theorem [1] is that cloaks may simultaneously comprise PM stiffness and 
anisotropic inertia, which provides a vastly richer potential set of material parameters, not 
limited to the model of eq. ([8|). 

Theorem [1] implies that the phase speed, wave velocity vector and polarization (not nor- 
malized) for plane waves with phase direction n are, from eq. ()38|) and (I39p . 

?;2 = n-V2n, c = v^^V'^n, q = S^^c. (44) 

The phase speed and wave velocity are independent of whether the cloak is an IC or the 
generalized PM-IC. These important wave properties are functions of the deformation only. 
They can be expressed in revealing forms using the deformation gradient as v = |F*n| and 
c = FN, where N = F*n/|F*n|. Note that the polarization q does in general depend upon 
the PM properties through the stress S. 

4.1.1 Continuity between the cloak and the acoustic fluid 

Continuity conditions at the cloak outer surface in the physical description follow in the 
same manner as (jlSp . The main difference is that the stress in the cloak is not isotropic, and 
therefore the condition that the shear tractions on the boundary vanish must be explicitly 
stated. The conditions for the pseudo-pressure which satisfies eq. (I36p are 

[nSp] =0, [n • p~^SVp] = on duj+. (45) 

These follow from eqs. (I33p and (I32p . 

4.1.2 Rays in the cloak are straight lines in the undeformed space 

Although Theorem [1] implies that the simple wave equation ([9]) in is exactly mapped to eq. 
(|36p in uj and hence all wave motion properties transform accordingly, including rays, it is 
instructive to deduce the ray transformation separately. We now demonstrate explicitly that 
rays in the cloak to, which are curves that minimize travel time, are just straight lines in 
Consider the straight line: X(t) = Xq + rN, where N is a unit vector in Q. The associated 
curve in oj is x(r) = x(Xo + rN). Differentiation yields dx/dr = FN, or dx/dr = V^s 
where the vector s is defined s = F^*N. Differentiating s(r), keeping in mind that N is fixed, 
gives 

^ = ^F*s = sV2f-*(^ • V)F-* = lsV2(VV-2)^, (46) 
dr dr Mr ^ 2 ^ Mr' ^ ^ 

where the compatibility identity OFf^^ /dxk = dFj^ /dxj has been used. We therefore deduce 
that straight lines in are mapped to solutions of the coupled ordinary differential equations 

i^.V^, ^^-is^lVVV (47, 
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But these are identically the ray equations in the cloak, see Appendix B. They are also the 
geodesic equations for the metric V~^. The ray equations conserve the quantity s- V^s which 
is equal to unity, reflecting the fact that s is the slowness vector, s = n/v, see eqs. (|44p and 
()B.4p . An illustration of rays inside the physical cloak is presented in Section [5j 

4.1.3 Relation to the Milton, Briane and Willis transformation 

Milton, Briane and Willis [10] examined how the elastodynamic equations transform under 
general curvilinear transformations. They showed, in particular, that if the deformation is 
harmonic then the constitutive relation ([6]) and momentum balance ([7]) for a compressible 
inviscid fluid with isotropic density transform into the equations for a pentamode material 
with anisotropic inertia, eqs. ()3ip and (I32p . respectively. The deformation is harmonic if 
V^x = 0, which realistically limits the transformation to the identity [10]. This would appear 
to indicate acoustic cloaking using the transformation method is impossible, in contradiction 
to the present result. In fact, as we show next, the MBW result is a special case of the more 
general theory embodied in Theorem [H one that corresponds to the choice S = J~^V^. 

The PM stiffness and inertia tensor found by Milton et al. [10] are 
and p = J~^V^ (their eqs. (2.12) and (2.13)). These are of the general form required by 
eq. (|43bp if we identify S as S = J~^V^. Does this satisfy the equilibrium condition (j43cp ? 
Using eq. ([5]), divS = div J^^FF* = J^^DivF* and this vanishes iff the deformation is 
harmonic. The MBW transformation therefore falls under the requirements of Theorem [1] for 
the specific choice of S = J~^V^ which satisfies the equilibrium equation (j43cp only if the 
deformation is harmonic. 

Having shown that the MBW transformation result is a special case of the present theory, 
it is clear that the transformation as considered here is different from theirs. Milton et al. 
|10j demand that all of the equations transform isomorphically, whereas the present theory 
requires only that the scalar acoustic wave equation is mapped to the scalar wave equation 
for the PM, see eqs. (j43ap . The mapping contains an arbitrary but divergence free tensor S 
which defines the particular but non-unique constitutive relation ([6]) and momentum balance 
([7]). Consider, for instance the displacement fields U(x) and u in and w, respectively. 
Under the transformation of [lO] xi(^x) u(x) = F~*U(x)(X) (eq. (2.2) of [10]). There is no 
analogous constraint in the present theory. In other words, we do not require an isomorphism 
between the equations for all of the field variables. Instead, the scalar wave equation for the 
acoustic pressure is isomorphic to the scalar equation for the pseudo-pressure of the PM. 

4.2 Cloaks with isotropic inertia 

Theorem [1] opens up a vast range of potential material properties. It means that there is no 
unique cloak associated with a given transformation O — > u; and its deformation gradient F. 
We now take advantage of this non-uniqueness to consider the possibility of isotropic inertia. 
Equation ()43bp indicates that the density is isotropic if S is proportional to V. Hence we 
deduce 

Lemma 4 A necessary and sufficient condition that the density is isotropic, p = pi, is that 
there is a scalar function /i(x) such that 

div(/iV) = 0, (48) 

in which case 

p = h^J, K = J, S = hY, (49) 
and the Laplacian is V^p = /iJV : V (^h^^ J^^W/p^ . 
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There is a general circumstance for which a solution can be found for h. It takes advantage 
of the second order differential equality 

V^p = (F*V) • F*Vp. (50) 

Although F is generally unsymmetric, F = F* in the special case that the deformation 
gradient is a pure stretch with no rotation (R = I). We therefore surmise 

Lemma 5 If the deformation gradient is a pure stretch (R = I and hence F coincides with 
then the density is isotropic, 

p = J-\ K = J, S = J"^V, (51) 

and the Laplacian becomes V^p = V : V ( VVp) . 

The infinite mass problem of the IC can be avoided if the material near the inner bound- 
ary duj- has integrable mass. This could be achieved, for instance, by requiring that the 
deformation near du- is symmetric (pure stretch) . Lemma [5] and the scaling arguments of 
Subsection 12.41 imply that the isotropic density scales as p =0(e'^~^), which is integrable as 
long as f3 < d + 1 {a> l/{d + 1)). 



4.3 Example: the rotationally symmetric cloak 

We again consider the deformation of eq. ()16p for the cloak = {x : < a < |x| < 6}, and 
assume the symmetric tensor S has the form S = w{r)(lr + 'y{r)l±). Differentiation yields 
divS = ^w' — (d — 1)(7 — l)]x, and the "equilibrium" condition ([43c]) is satisfied if tt;(r) and 
7(r) are related by w' = {d — 1)(7 — 1). It is convenient to introduce a new function g{r) 
such that 7 = rg'/g and w = {g/r)'^~^, which automatically makes divS = 0. The cloak 
parameters therefore have general rotationally symmetric form 



/ „2 

^ ^ (52) 



The functions / and g are independent of one another, and together define a 2-degree 
of freedom class of PM-IC cloaks. The general solution has both anisotropic stiffness and 
anisotropic inertia. The previous example of the pure IC corresponds to the special case of 
g = r, for which eq. (f52|) gives S = I and K and p agree with eq. (fT8]l . 

The form of the stress S indicates the PM-IC has transversely isotropic (TI) symmetry. 
This is a special case of the orthotropic PM considered earlier. A normal TI solid with axis 
of symmetry in the X3— direction has five independent elastic moduli: Cn, C33, C12, C13 and 
C44. The last is a shear modulus, the other shear modulus is Cge = (C'n — Ci2)/2. We set 
all shear moduli to zero - implying C44 = and C12 = Cn, and the remaining independent 
moduli Cii, C33 and C13 satisfy C11C33 — Cf^ = 0. The PM therefore has two independent 
elastic moduli. Let C33 — > Kr{r), Cu K±(r) and C13 — > y/KrK±, then the fourth order 
elasticity tensor defined by (152]) is 



C = KS(g)S = {yjKrlr + VK_^I±) ® (V^Ir- + V^I±)> (53) 



where the stiffnesses K^, K±, and the principal values of the inertia tensor given by eq. (|19p. 
are 

Kr = Jj{^) , I^± = JJ7{^) , Pr = f{^) , P± = ^{^) . (54) 



The phase speeds and c± in the principal directions are again given by eq. (j2ip . This 
might seem amazing at first sight, but recall that it is predicted from the general theory. 
That is, the phase and wave velocity speeds are independent of how we interpret the cloak 
material, as an IC or the more general PM-IC. In the present example, it means that the 
phase and wave velocity are independent of g. 
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Table 2: Behaviour of quantities near the vanishing point r = a for the scahng / oc as 
^ = r — a I with isotropic inertia. 



dim 




c± 


Kr 


K± 


P 


2 












3 








^a— 1 





4.3.1 Pure PM cloak with isotropic density 

The inertia is isotropic when pr = p±, which occurs if g{r) = f{r). In that case p = pi, and 
eq. reduces to 

P = n^Y-\ K,- = j{^r\ K^ = f'{^Y-\ (55) 

We observe that the parameters of eq. (j55p are obtained from the IC parameters in eqs. 
([18]) and (fT9]) under the substitutions {K, p^, p±} — > {1/p, l/i^r, ^/K±}. Thus, the universal 
relation analogous to eq. ([20]) is now 

KrKi-^=p''-\ (56) 

and by analogy with eq. ([22]) the three original material parameters can be expressed using 
the phases speeds only, as 

p = c-^(H'i-^\ Kr = CrC-^'^-^\ K^ = c-^d^-\ (57) 

In summary, there is a one-to-one correspondence between the two sets of three material 
parameters for the limiting cases of the pure IC on the one hand, and the pure PM cloak on 
the other. Of course, as discussed before, the density and stiffness cannot be simultaneously 
isotropic. The PM-IC cloak with material properties (|52p includes both limiting cases when 
g = r and g = f, respectively. 

Table 2 summarizes the scaling of the physical quantities for isotropic inertia, similar to 
the scalings in Table 1 for the pure IC. Note that the wave speeds c,., c± and the intermediate 
(C13) modulus \/KrK± have limiting behaviour which is independent of the dimensionality, 
while the density p and the moduli Kr and K± depend upon whether the cloak is in 2D or 
3D. 



5 Further examples 

5.1 A non-radially symmetric cloak with finite mass 

The examples considered above are rotationally symmetric and rather special in that they 
can be made using uniformly pure IC, or pure PM, or hybrid PM-IC. The pure IC model is 
always achievable as Lemma [1] showed, but it suffers from the infinite mass catastrophe. The 
pure PM model requires that Lemma [U hold at all points, which is not realistic. However, 
we can always obtain a cloak comprising partly pure PM by requiring the deformation to be 
locally a pure stretch (Lemma [5]). In particular, by constraining the deformation near the 
inner surface duj- in this manner, the density can be made both isotropic and integrable. We 
now demonstrate this for a non-rotationally symmetric cloak. 

For A eSym"*" and h{Q, h'{() > for ^ £ [0, 1], consider the deformation 

X = C"^/i(C)AX, C = (X-AX)^/l (58) 
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This generalizes the deformation of eq. ()16p (A = I) and has the important property that 
the deformation gradient is symmetric, 



F = ^A + -(-)' (AX)® (AX). (59) 

The inner surface is an elhpse (2D) or eUipsoid (3D), 

5cj„ = {x : X • A-^x = /i^(0)}. (60) 
The mapping X <-> x must be the identity on the outer surface of the cloak (9u;+ = 90. 



This eliminates the transformation (|58p as a possible deformation in the vicinity of dto-^- but 
it does not rule it out elsewhere. In particular, it can be used on the inner surface duj- and 
for a finite surrounding volume. Then it could be patched to a different mapping closer to 
the outer boundary of the cloak, one which reduces to the identity on duj^. For instance, 

X = C'hiOA'^X, (61) 

where J^(x) = 1 for all x between duj- and some surface C, beyond which v decreases smoothly 
to zero as x approaches dto^, which is assumed to be a level surface of i.e. an ellipsoid or 
an ellipse. We assume that C = 1 on the outer surface, so that 

dLO+ = {x: x-Ax = h^{l)}. (62) 

Let C be the level surface C = Co foi' constant Co ^ (O;!)- The surface separating the pure 
PM inner region from the PM-IC outer part of the cloak is therefore 

C = {x: X- A-^x = h^iCo)}. (63) 

Based on Lemma [5] the inner part of the cloak between duj^ and C can be constructed from 
pure PM material with isotropic density p = = {\A\h')~^ (^(^/h^ ^. The remaining part 
of the cloak is PM-IC, and the mass of the entire cloak will be finite. 

For instance, in Figure E) h{C) = ^(1 + C) for C e [0,1]; = 1 for C G [0, |) and 
u = 4(1 — ^) for C G [|, 1], and the principal values of A are 0.6 and 1.0. The Figure shows 
each ray following a continuous path through the cloak, with collinear incident and emergent 
ray paths. There is a unique ray separating the rays traversing the cloak in opposite senses, 
and which defines a "stagnation point" at the cloak inner surface. The separation ray is the 
one that would intersect the singular point in the undeformed space, O in Figure [H This 
is the origin in Figure [2] and since the rays are incident horizontally, the separation ray is 
defined by rE2 = outside the cloak, and it intersects doj- at x = ±^Aei/ y/ei ■ Aei. The 
wavefront in effect splits or tears apart at the incident intersect and it reforms at the emergent 
intersect. The time delay between these two events is infinitesimal since the tearing/rejoining 
is associated with the instant at which the wavefront would traverse O in the undeformed 
space. A time-lapse movie illustrating this may be seen here (20 seconds long). Another 
movie showing the ray paths for different directions of incidence can be found at the same 
URL. 

5.2 Scattering from near-cloaks 

A near-cloak or almost perfect cloak is defined here as one with inner surface dcu^ that does 
not correspond to the single point X = O. We illustrate the issue using the radially symmetric 
deformation (I16p with /(a) small but non-zero, and assuming time harmonic motion, with 
the factor e"*'^* understood but omitted. Since the inner surface is not the image of a point, 
it is necessary to prescribe a boundary condition on the interior surface, which we take as 
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1.5 





Figure 2: Ray paths through a non-radially symmetric cloak. The sohd curves are the inner 
and outer surface of the cloak. The dashed curve delineates the inner region in which the 
deformation gradient is everywhere symmetric and the cloak is pure PM with finite isotropic 
mass. Two movies of the rays and wavefronts in this cloak may be viewed Iherel 



zero pressure on r = a. The specific nature of the boundary condition should be irrelevant 
as /(a) shrinks to zero. 

As before, the cloak occupies uj = {x : a < r < b}, but now /(a) > 0. The total response 
for plane wave incidence is p(x) = pQ^^^Ri"^)^"^^ -|-psc(x), where po is a constant and 

°o fi"(2-(5„o)Jn(fci?(a))^^C^^^cosne, ( f(^\ „ < r <^ h 

^0 zn2n + l)j„(fci?(a))J|^ygf:^p„(cose), l^-, 6 < r < oo, 

^. h'„'(kR{a)) 

in 2D and 3D, respectively. A near-cloak can be defined in many ways: for instance, a power 
law /(r) = b{^y with 0<6 < a is considered in [22]. Here we assume a linear near-cloak 
mapping similar to the one examined by Kohn et al. [20j . 



= , (64) 

where < 6 < a. Hence, f^^\a) = 5 and the radius at which the mapping is zero, r = 
a — 6(b — a)/{h — 5), defines the size of a smaller but perfect cloak. 

Some representative results are shown in Figure [3l which illustrates clearly a disparity 
between cylindrical and spherical cloaking, even when the physical optics cross-sections are 
identical. Thus, for /(a) = O.Ola the 3D cross-section is negligible (Fig. [3ld) but for 2D the 
cross-section is two orders of magnitude larger (Fig. [3lc). Ruan et al. [23] found that the 
perfect cylindrical EM cloak is sensitive to perturbation. This sensitivity is evident from the 
present analysis through the dependence on the length 6 which measures the departure from 
perfect cloaking (5 = 0). 

The ineffectiveness of the same cloak in 2D as compared with 3D can be understood 
in terms of the scattering cross-section. The leading order far-field is of the form p = 
Poe^^^ + Po5(x)r~('^~-^)/^(i27r/A;)(^"'^)/^e*'^''. The optical theorem implies that the total scat- 
tering cross-section, and hence the total energy scattered, is determined by the forward 



15 



(a) 2D, S = 0.3, E = 1.48 



(b) 3D, 5 = 0.3, E = 0.79 




(c) 2D, S = 0.01, E = 0.12 



(d) 3D, S = 0.01, E = 1 X 10" 



Figure 3: A plane wave is incident from the left with frequency = 10 on the cloak defined by 
eq. (|64p with a = 1, b = 2^/^'^~^\ The outer cloak radius b is chosen so that the geometrical 
cross-section of the cloak is twice the geometrical cross-section of the cloaked region in both 
2D and 3D. The circular core in the plots is the cloaked region of radius a. The virtual inner 
radius /(a) = 5 is 0.3 or 0.01, and S is the total scattering cross-section. 
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scattering amplitude: S = 47rA; ^l'nig{ez). Thus, 



I V (2 - 5„o)Re4#^^, 2D, 



# g (2n + l)Re^#^, 2,D. 

The cross-section is dominated in the smah kf[a) hmit by the n = term, with leading order 
approximations 



^|lnA:/(a)|-^ + ..., 2D, 
47r/2(a) + ..., 2,D. 



(66) 



This explains the greater efficacy in 3D, and suggests that all things being equal, cylindrical 
cloaking is more difficult to achieve than its spherical counterpart. 



6 Discussion and conclusion 

Starting from the idea of an acoustic cloak defined by a finite deformation we have shown 
that the acoustic wave equation in the undeformed region is mapped into a variety of possible 
equations in the physical cloak. Theorem [1] implies that the general form of the wave equation 
in the cloak is 

KS : V(p"^SVp) - p = 0, (67) 

where the stress-like symmetric tensor S is divergence free and the inertia tensor is p = 
JSV~^S. The non-unique nature of S for a given fixed deformation opens many possibilities 
for interpreting the cloak in terms of material properties. 

If S is constant (S = I with no loss in generality) then the cloak material corresponds to 
an acoustic fiuid with pressure p defined by a single bulk modulus but with a mass density 
p that is anisotropic, which we call the inertial cloak (IC). The IC model is mathematically 
consistent but physically impossible because it requires a cloak of infinite total mass. There 
appears to be no way to avoid this if one restricts the cloak material properties to the IC 
model. If one is willing to use an imperfect cloak with finite mass, and is concerned with fixed 
frequency waves, then the scattering examples show that significant cloaking can be obtained 
by shrinking the effective visible radius to be subwavelength. The 2D and 3D responses for 
imperfect cloaking are quite distinct, with far better results found in 3D. 

A cloak of finite mass is achievable by allowing S to be spatially varying and divergence 
free. The general material associated with eq. (j67p . called PM-IC, has both anisotropic inertia 
and anisotropic elastic properties. The elastic stiffness tensor has the form of a pentamode 
material (PM) characterized by the symmetric tensor S and a single modulus K. Under 
certain circumstances, characterized in Lemmas H] and O the density becomes isotropic and 
the material is pure PM. More importantly, the total mass can be made finite. 

The finite mass problem arises from how we interpret the cloak material in the neigh- 
bourhood of its inner surface. It is therefore not necessary to totally abandon the pure IC 
model, but it does mean that the alternative PM-IC is required at the inner surface. Prom 
the examples considered here it appears that one can always use a pure PM model near the 
inner cloak surface, and thereby achieve finite mass. One method is to force the deformation 
near the inner surface to be a pure stretch, then Lemma [5] implies that the density is locally 
p = 1/det F. The total mass remains finite as long as p is locally integrable, which is easily 
achieved. 

The theory and simulations of PM-IC and PM materials presented here illustrate the 
wealth of possible material properties that are opened up through the general PM-IC model 
of acoustic cloaking. Physical implementation is in principle feasible: for instance, anisotropic 
inertia can be achieved by microlayers of inviscid acoustic fiuid [Tl] , while the microstructure 
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required for pentamode materials has been described [18]. It remains to combine these known 
methods to obtain practical PM-IC materials. 
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Appendix 

A Proof of Theorem [1] 

A weak but instructive form of the identity (|41|) is proved first. Consider the possible identity 

V^p = cA: V(B~2divc-^Ap), (A.l) 

where A, B GSym are non-singular and c is a scalar. Let us examine under what circum- 
stances this identity holds. Let q be an arbitrary test function, and consider the integral 

1= f dygcA : V(B~2divc"^Ap). (A.2) 
Jn 

Substituting dV = dv, integrating by parts and ignoring surface contributions, yields 

I = - / d t; ( div J'^cAq) ■ (B'^ div c'^Ap) . (A.3) 

In order to guarantee the integral is self-adjoint, that is, symmetric in both p and we 
demand c = J^/^. The self-adjoint property is made evident by writing / as 

I = - I dv (B-^ div J-^'^Aq) ■ (B-i div .r^'^Ap) . (A.4) 

If (|A.ip is to be vahd, then 



/=- / dVVxq-^xV- (A.5) 
Comparing these integrals and once again using dv = J dV , implies 

( Ji/^g-i j-i/2Ag) . ( j1/2b-i div ,r^/^Ap) = Vxq ■ ^xP- (A.6) 
The only way that these can agree for arbitrary p and q is if 

divJ~i/2A = 0, (A.7) 

in which case ()A.6p becomes 

(B-iAF-*Vxg) • {B-^AF-'Vxp) = V xQ ■ ^ xP- (A.8) 

Using = FF*, it is clear that (fXS]) can only be satisfied if 

B-2 = A-1v2a-^ (A.9) 

A weak form of Theorem [1] follows by substituting A = J^/^P. Based upon this it is a 
straightforward exercise to see that the identity ()4ip can be derived directly by brute force 
differentiation of the right hand side, taking into account the constraint (j42p and Lemma [TJ 
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B Ray equations in an acoustic cloak 

Consider a WKB type of solution for the displacement: u(x, t) = U(x, t)e'^'"^^^'^\ The leading 
order equation for the phase cp and amplitude U is (see eqs. (p3|) - (j38|) ) 

[K(SV(/))®(SV0)-0V]U = O. (B.l) 

The inner product of eq. (jB.ip with U may be written H^H^ = where 

H± = ^±K^/^\ci- p-^/^SV(l)\, (B.2) 

and q = p^^'^XJ /\p^^'^\J\. We focus on the characteristic H^((j>,\7(p,x) = 0. The Hamilton- 
Jacobi equations for this "Hamiltonian" yield the ray equations, 

dt _ dH+ dx _ dH+ dVcp _ _dH+ d(/> _ _dH+ 
d^~~df' d7~9V^' ~ 5^' d7 ~ dT' ^ 

where r is the time-like ray parameter. Since dH^/d(j) = 1, r may be replaced by t as the 
natural ray parameter, while dH^ /dt = implies that (j) is constant along a ray. We choose 
(p = —1 for convenience. Define the slowness vector s(t) = along the ray x = x(r). The 
vector equation ()B.1|) then implies that q = i^^/^p~ ^/^Ss, and since q is by definition a unit 
vector, using eq. (|43bp we deduce that the slowness satisfies 

s • V^s = 1. (B.4) 

This is simply the ellipsoidal slowness surface mentioned in Section HI Finally, the evolution 
equations along the ray can be expressed as a closed system for x(t),s(t) by using eq. (|B.3p 

and noting that = (s • V^s)^''^ — 1, 

^=V^s, ^ = -(Vs) . (VV)s. (B.5) 
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